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Abstract. We review various prescriptions which have been proposed for the turbulent transport of matter and angular momen- 
tum in differentially rotating stellar radiation zones. A new prescription is presented for the horizontal transport associated with 
the anisotropic shear turbulence which is produced by the differential rotation in latitude; this '/?- viscosity' is drawn from torque 
measurements in the classical Couette-Taylor experiment (Richard & Zahn 1999). Its implementation in a stellar evolution code 
leads to enhanced mixing, as illustrated by models of a rotating main-sequence star of 1.5 solar mass. 
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1. Introduction 

In standard models of stellar interiors, radiation zones, which 
are convectively stable, are postulated to be without motion 
other than rotation, except in the vicinity of convection zones, 
where often convective overshoot is accounted for in a crude 
way. But it is well known, since Eddington (1925) and Vogt 
(1925), that in rotating stars the centrifugal force breaks the 
radiative equilibrium, and thus causes a slow large-scale cir- 
culation, which transports matter and angular momentum. A 
similar circulation is induced locally in the vicinity of a differ- 
entially rotating convection zone, in the so-called tachocline. 
As a result, radiation zones rotate non-uniformly, and therefore 
they are liable to various instabilities which may generate tur- 
bulence. 

This turbulence contributes to the transport of matter and 
angular momentum, and to model these in a stellar evolution 
code, one needs reliable prescriptions for the turbulent trans- 
port. To derive such prescriptions remains today one of the 
most challenging problems, not only in stellar physics, but also 
in oceanography and in atmospheric sciences. The reason is 
that turbulence is still poorly understood, particularly when it 
departs from ideal conditions: homogeneity, isotropy, station- 
arity. "Turbulence is an order of magnitude subject; there is no 
rigorous theory of turbulence" warned M. Hoffert at the Aspen 
Global Change Institute in 1997. And he concluded: "The way 
turbulence is represented in current (global circulation) mod- 
els is not realistic and this has important ramifications for the 
climate sensitivity predicted by these models." 
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We could use the same words to describe the situation when 
modeling stellar interiors, replacing "climate sensitivity" by 
"evolutionary sequences". Therefore it takes some courage to 
- nevertheless - engage in this task, as we shall attempt in this 
article. We shall restrict our scope on turbulence produced by 
shear instabilities, which are likely to be the most powerful, at 
least in a non-magnetic star. 

Why do we believe that stellar radiation zones are turbu- 
lent? Some velocity profiles are known to be linearly unstable, 
as soon as the Reynolds number characterizing the flow ex- 
ceeds some critical value. (We recall that this number involves 
the typical size L and and the typical velocity contrast V of the 
flow: Re = VL/v, where v is the kinematic viscosity; it mea- 
sures the relative strength of advection and viscous transport.) 
But flows that are linearly stable are found turbulent in labora- 
tory experiments, again above some critical Reynolds number 
Re c , due to non-linear processes which are still not completely 
elucidated (see for instance Richard & Zahn 1999). Owing to 
their huge size and their low viscosity, stellar radiation zones 
have a Reynolds number which easily exceeds the critical val- 
ues measured in the laboratory, Re » Re c , and therefore they 
should be turbulent, according to this simple criterion. 

What seems more powerful in preventing a shearing flow 
from becoming unstable is the buoyancy force in a stable strati- 
fication, or the Coriolis force for fast enough rotation. We there- 
fore expect these forces to play an important role in controlling 
the turbulence in stellar radiation zones; in particular, the turbu- 
lent motions will probably be anisotropic, with stronger trans- 
port in the horizontal than in the vertical direction. 

Several papers have been devoted to the description of such 
shear turbulence, and to the derivation of prescriptions for the 



2 



S. Mathis et al.: On shear-induced turbulence in rotating stars 



turbulent transport; to quote the most recent ones we refer to 
Zahn (1992), Maeder (1995, 1997), Talon and Zahn (1997), 
Schatzman et al. (2000), Maeder (2003). Here we shall present 
an update of the subject, and recall the main results. We shall 
also introduce a new prescription for the horizontal transport, 
which is drawn from laboratory experiments, and discuss its 
impact on a stellar model, namely that of an evolving 1 .5 M 
star. 

2. Turbulence induced by the vertical shear 

It is well established that a stable entropy stratification acts to 
inhibit the instability caused by a vertical shear V(z), where V 
is the amplitude of the horizontal velocity and z the vertical 
coordinate. 

In the adiabatic case, i.e. in the absence of thermal leakage 
and viscous dissipation, a necessary condition for shear insta- 
bility is expressed by the Richardson criterion: 



(1) 



It states that the shear rate dV/dz has to overcome the stabi- 
lizing effect of the buoyancy force, which is measured by the 
Brunt-Vaisala frequency A^: 



N 2 = N 2 + Nl = j- (V ad 
tip 



■ V) + — v„ 

Hp * 



(2) 



We employ the usual notations for temperature T, density p, 
pressure P, molecular weight p, gravity g, pressure scale height 
Hp, and for the logarithmic gradients V = dlnT/dlnP and 

= dln/v/dlnP; V a d = {d\nT /<91nP) ac | is the adiabatic tem- 
perature gradient. The coefficients 5 - -(d In p/ din T)ph and 
(f = (d In p/d In p)pj are drawn from the equation of state; they 
are unity for perfect gas. The critical Richardson number Ri c 
is of order unity for typical flow profiles; in the following we 
shall take Ri c = 1 /4. 

However, when the turbulent eddies exchange heat with 
their environment, the stabilizing effect of buoyancy is re- 
duced, as was shown by Townsend (1958). In an optically thick 
medium, this reduction is determined by the Peclet number 
characterizing the turbulence: 



Pe = 



vt 



where v and t are the velocity and the size of the largest turbu- 
lent elements, and K the thermal diffusivity: 



K = 



16 err 3 
3 Cpp 2 K 



cr is the Stefan constant, k the Rosseland mean opacity and Cp 
the heat capacity at constant pressure. 

When vt <sc K thermal diffusion proceeds faster than ad- 
vection and, if there is no composition gradient, the Richardson 
criterion for instability takes the modified form (Zahn 1974 1: 



vt , (dvy 



(3) 



Identifying the largest vt satisfying this inequality with the tur- 
bulent diffusivity D,, one obtains 



D, 



K RUdV 
Nl \ dz 



(4) 



This result is consistent with the linear stability analysis of a 
hyperbolic-tangent profile imbedded in a stable stratification, 
including thermal diffusion, which were performed by Dudis 
( 1974) and more recently by Lignieres et al. (1999). 

In presence of a vertical composition gradient, only the 
thermal component of the buoyancy force is weakened through 
radiative diffusion, and the instability criterion would read then 
(Zahn 1974; Maeder 1995) 



vt , , dV 
:Nl + N 2 < Ri, 



K 



dz 



(5) 



Meynet and Maeder (1996) pointed out that this condition was 
too severe, and that it would not allow for the mixing which is 
observed in massive stars. 

But the //-component of the buoyancy can also be reduced 
by the turbulence. This is true in particular when this turbu- 
lence is anisotropic, with much stronger transport in the hor- 
izontal than in the vertical direction, as expected here due to 
the stable stratification. Then v and t characterize the largest 
eddies which are still isotropic, and the horizontal component 
of the turbulent diffusivity is much larger than the vertical one: 
D/, » vt. It was argued by Talon and Zahn (1997) that such tur- 
bulence erodes the fluctuations of chemical composition which 
provide the restoring force, and that the Richardson criterion 
then takes the form 



vt 



vt 



(K + D h ) 



-m + — Nf, < RL — 



D h 



dV 



dz 



(6) 



The rotation rate Q. does not appear in this inequality because 
the buoyancy force is stronger than the Coriolis force, N ?s> Q.. 
From (|6} it is straightforward to derive the vertical component 
D v = vt of the turbulent diffusivity: 



D v 



Ri c 



N 2 T l{K + D h ) + NllD h 



dr / 



(7) 



we have replaced the shearing rate by its expression in spherical 
geometry: dV/dz — » rsin#dQ/d#, and performed a spherical 
average. The critical Richardson number has been renormal- 
ized; its canonical value is now Ri c = 1/6. A similar expression 
is being used by Maeder and Meynet (2001). 

3. Turbulence induced by the horizontal shear 

Differential rotation in latitude, Q(0), is known to be linearly 
unstable when one of the following conditions is met. Either 
the specific angular momentum (r sin 0) 2 Q has to decrease from 
pole to equator, on a level surface, according to the H0iland- 
Solberg criterion, a situation which is highly improbable in 
stars. Or the vorticity profile <9(sin 2 6>Q)/ sin 986 must present 
a maximum in 0, a criterion akin to the famous inflexion 
point theorem (Rayleigh 1880), which has been transposed by 
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Watson ( 1 98 1 ) to spherical geometry. This requirement is prob- 
ably not satisfied in stars, except perhaps in a tachocline (see 
Garaud2001). 

But in stars the Reynolds number characterizing such dif- 
ferential rotation is easily so large that it is tempting again to 
conclude that even a slight shear, although linearly stable, is 
prone to non-linear instabilities. However some caution must 
be exerted here, since the Coriolis force has a stabilizing effect 
on shear instabilities, and the role of curvature remains to be 
fully elucidated. Although work is in progress on that subject, 
most of what is known so far is based on laboratory experi- 
ments. There turbulence tends to suppress the cause of the in- 
stability, namely the differential rotation when the turbulence is 
generated by the shear (Richard & Zahn 1999). We shall there- 
fore assume that the same occurs in stellar radiation zones, and 
that the main role of the turbulent viscosity v/, is to reduce the 
horizontal shear produced by the advection of angular momen- 
tum through the meridional circulation. This is described by the 
advection/diffusion equation governing the differential rotation 
in latitude: 



p^-(pr 2 Q 2 ) - 2pQr[2V 2 - aU 2 ] 



(8) 



10pv/,Q 2 , 

where the meridional circulation is given by 

u r (r, 9) = U 2 (r)P 2 (cos 9), u e {r, 9) = V 2 (r)dP 2 /d9, (9) 

W ith^ 2 (r) = ±±(p r*U 2 ), a(r)= 1 -^*, (10) 
opr ar v ; 2 dlnr 

and where we have expanded the angular velocity in Q(r, 9) = 
Q(r) + Q 2 (r)[P 2 (cos ff)+ 1/5] 1 , the mean angular velocity being 
defined as 

_ r Q(r, 9) sin 3 9d9 

ft(r)= „ ■ (11) 

J Q sin 3 6»d6» 

In the stationary limit, advection and diffusion balance each 
other: 



Qr[2V 2 - aU 2 ] = 5v h D. 2 . 



(12) 



For lack of a better prescription, it was proposed in Zahn 
(1992) that v/, be determined by the condition that it enforces a 
small, postulated amount of differential rotation Q 2 /Q. Then 

v h = -r[2V 2 -aU 2 ], (13) 

Ch 

where Cy, < 1 is a parameter of order 1, with a similar expres- 
sion for the turbulent diffusivity £)/,. 

The horizontal diffusivity due to that shear turbulence in- 
terfere with the transport of chemicals through the meridional 
circulation, as was described by Chaboyer and Zahn (1992). 
When the turbulence is sufficiently anisotropic, the combina- 
tion of both processes (meridional advection + horizontal dif- 
fusion) results in a vertical diffusion whose coefficient is given 
by 



£>eff = 



1 (rU) 2 
30 D[, 



(14) 



1 The constant 1/5 is missing in Zahn (1992) and Maeder & Zahn 
(1998); the correct expression is derived in Mathis & Zahn (2004), as 



and which competes with that, D v » v v , caused by the vertical 
shear. 

As Maeder (2003) has shown very convincingly, the crude 
prescription ( II 31 given above allows for horizontal fluctuations 
of the molecular weight which are too large and tend to stop the 
meridional circulation. It also leads to a vertical diffusivity D v , 
according to 0, which is insufficient to account for the mixing 
observed in massive stars. Moreover, it does not guarantee that 
Dh » D v , as is assumed in deriving the effective diffusivity 
Deff. 

For these reasons, Maeder proposes another prescription 
for this horizontal shear viscosity, which is based on the bal- 
ance between the rate at which kinetic energy is channeled into 
differential rotation by the horizontal component of the merid- 
ional circulation, and the rate of viscous dissipation. He as- 
sumes that the dynamical time involved in the first process can 
be approximated by the time needed by a lagrangian particle 
drifting from equator to pole to be wound up by (say) 27r. This 
yields a horizontal viscosity 

v h = Ar(rQV 2 [2V 2 - aU 2 ]) u \ (15) 

with A of order 1/10. 

Here we take another approach, and borrow a prescrip- 
tion for the turbulent transport from laboratory experiments. 
Cylindrical Couette-Taylor flow has been found turbulent in the 
laboratory when the rotation rate O increases with the distance 
s from the axis, whereas the linear instability analysis predicts 
that the flow should be stable. Experiments performed with the 
inner cylinder at rest, and varying the size As of the gap be- 
tween cylinders, have shown that when the relative gap As/s is 
larger than about 1 /20, turbulence sets in whenever the gradi- 
ent Reynolds number s 3 (d£2/ds)/vexceeds some critical value. 
A turbulent viscosity can be drawn from the torque measure- 
ments, and it is found to scale as 
dC2 



ds 



(16) 



with p * 1.5 10" 5 (Richard & Zahn 1999). This '/3-viscosity' 
has been applied to accretion disks by Hure et al. (2001), where 
it yields somewhat different results compared to the classical 
a-viscosity introduced by Shakura & Sunyaev (1973). 

This prescription has been established in the case of maxi- 
mum differential rotation; it remains to be verified whether its 
validity extends to milder shear rates, which are typical of stel- 
lar interiors. Transposing d!6i to spherical geometry, with the 
shear in latitude only, 



v, =p(r sin 6>) 



dQ. 



rd9 



:p> 2 |Q 2 |sin 3 i 



dPi 



de 



(17) 



and averaging over the sphere, we obtain a turbulent diffusivity 
1 



-p r l \n 2 



(18) 



Inserting this value in (I12> . we get the following expression for 
the turbulent viscosity 



Vh 



1/2 



(r £ D.) 1/2 [r\2V 2 - aU 2 \] 



1/2 



(19) 



well as eq. JHJ. 



again in the stationary limit. The same value will be used for 
the turbulent diffusivity D/,. 
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this work Maeder (2003) 

Zahn (1992 



M = 1.5 M s , Z=0.02, t = 368 Myrs 



this work Maeder (2003) 

Zahn (1992) 



M - 1.5 M s: Z-0.02, I - 1.36B Gyre 



0.4 0.6 
M/M. 




Fig. 1. Horizontal component D/, of the turbulent diffusivity 
generated by the horizontal shear, in a 1.5 M star rotating ini- 
tially with an equatorial velocity of 110 km s _1 , at two evolu- 
tionary stages (0.368 and 1.368 Gyr). Comparison of the old 
prescription (eq. Zahn 1992) with that of eg. 1151 (Maeder 
2003) and that of the present paper (eq.U9ll. 



Fig. 2. Same comparison for the effective diffusivity D e ff in the 
vertical direction, which results from the erosion of the advec- 
tive transport through the horizontal diffusivity Di, displayed 
in fig. 1. Here D e ff has been evaluated using the old prescrip- 
tion for D h (Zahn 1992), that of Maeder (2003) and that of the 
present paper (eq.fT^l. 



4. Comparing different prescriptions for the 
horizontal turbulent diffusivity 

To compare the different diffusivity prescriptions presented 
above, it is necessary to implement them in a stellar structure 
code, since they depend on quantities, such as the molecular 
weight gradient and the differential rotation, whose evolution 
in time is governed by these diffusivities. For this purpose we 
ran an evolutionary sequence of a 1.5 M star, of Population I 
composition (X=0.705, Y=0.275, Z=0.020), using the code 
STAREVOL. 

The transport of chemicals and angular momentum was 
treated according to Zahn (1992) and Maeder & Zahn (1998). 
The fourth-order advection/diffusion equation governing the 
evolution of angular momentum in the radial direction is split 
into a system of four plus one first-order equations of order one 
and it is solved through a Newton-Raphson method. Convective 
zones are supposed to rotate uniformly. We refer to Palacios et 
al. (2003) for more details concerning the physics used to com- 
pute of these models. 

The initial equatorial velocity was taken 110 km s _1 , 
and the star was submitted to magnetic spindown following 
Kawaler (1988) so that its equatorial velocity at the age of the 
Hyades ranges between 90 and 100 km s _1 , as expected from 
vsiru measurements in this cluster (Gaige 1993). The results 
are displayed in figs. 1-3, for two ages: 0.368 and 1.368 Gyr. 
In fig. 1 we compare the different prescriptions for the horizon- 
tal diffusivity £)/,, which is generated by the latitudinal shear. 
Note that the more recent prescriptions, that of Maeder (eq. H5t 
and that of the present paper (eq. [l9j, yield values which ex- 
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Lit = 1.5 M 3 , Z-0.02. L - 368 Myrs 




M = 1.5 M„, Z = 0.02, t = 1.36B Gyrs 



0.2 0.4 0.6 0.8 1 

M/M. 



Fig. 3. Same as in figs. 1 and 2 for the vertical component D v 
generated by the vertical shear. D v depends on D/ t when there is 
a vertical composition gradient, according to eq. [7J it has been 
evaluated here using the old prescription for D/ t (Zahn 1992), 
that of Maeder (2003) and that of the present paper (eq.U9>. 

ceed those of the old prescription (eq. 1 1 31 by 3 to 4 orders of 
magnitude. 

This enhancement of D/ t has two effects. First, it maintains 
the differential rotation in latitude at very low level; this is re- 
quired by the assumption of 'shellular rotation' made in cal- 
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culating the meridional circulation (Zahn 1992), namely that to 
first approximation the angular velocity may be considered as a 
function of r only. The second effect is to reduce the latitudinal 
fluctuations of chemical composition, which otherwise would 
tend to suppress the meridional circulation, as was pointed out 
already by Mestel (1953). The consequence is an increase of 
the circulation speed U2, which compensates that of Dh, and 
this explains why the effective diffusivity D e ff is less affected 
as one would expect by the increase of Dh, as can be seen in 
fig. 2. 

There are also qualitative changes in the profiles. With the 
Dh -prescription introduced here, the radial component U2 of 
the meridional circulation keeps everywhere the same (nega- 
tive) sign, and therefore D e ff never vanishes, as it does with 
the Maeder (2003) prescription, for which the circulation splits 
into three superposed cells. 

Finally, this larger value of D/, has the effect of increasing 
the turbulent diffusivity D v in the regions where it is reduced 
by a molecular weight gradient, according to @. This is illus- 
trated in fig. 3, and one can check that D v is everywhere four 
orders of magnitude smaller than Dh, which is consistent with 
the anisotropy assumption made in deriving D e ff. 



5. Conclusions 

We have presented here a new prescription, the viscosity', 
for the horizontal component Dh of the turbulent diffusiv- 
ity due the differential rotation in latitude, which is drawn 
from Couette-Taylor experiments. According to this prescrip- 
tion (eq. I19l l. Dh exceeds by 3 or 4 orders of magnitude that 
which was initially suggested (eq.H3> and which has been used 
until now. Remarkably enough, the prescription established by 
Maeder (2003) yields a value which is very similar, although it 
is based on quite different arguments. We have analyzed the ef- 
fect of this dramatic change in D/, on the structure of a 1.5 M 
star: it increases the vertical transport, and therefore enhances 
the mixing, particularly as the star evolves and a molecular 
weight gradient is building up. 

The main incentive of trying to improve this prescription 
came from the evolutionary calculations performed by Maeder 
and Meynet, who found that the turbulent transport, as it was 
modeled initially, was insufficient to account for the observed 
properties of massive stars. 

This new prescription may provide a better agreement with 
the observations, but we can hardly consider it as the final an- 
swer, given the uncertainties involved in applying to stellar in- 
teriors a result obtained in the laboratory, in cylindrical geom- 
etry and for extreme differential rotation. However we believe 
that this is the best we can offer today, until numerical sim- 
ulations allow us to model that type of turbulence in a more 
realistic way. 
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